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ABSTRACT 

We describe a theoretical framework to compute the cluster gas distribution in 
hydrostatic equilibrium embedded in a class of spherical dark matter halo potentials. 
Unlike the conventional isothermal /3-model, the present method provides a physical 
basis to directly probe the shape of dark matter halo from the observed X-ray surface 
brightness and temperature profiles of clusters of galaxies. Specifically, we examine 
the extent to which the resulting gas density and X-ray surface brightness profiles are 
sensitive to the inner slope of the dark matter halo density and other more realistic 
effects including the self-gravity of the gas and the polytropic equation of state. We 
also discuss a practical strategy to apply the present methodology to the actual cluster 
profiles from future X-ray observations. 



Subject headings: cosmology: theory - dark matter - galaxies: clusters: general — 
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1. Introduction 

The gas density profiles of X-ray clusters of galaxies are known to be approximated well by 
the empirical formula, the isothermal /3-model: 

= [l + (^/'°)2]3/3/2- (1) 

Theoretically this is consistent with the observed indication that luminous member galaxies obey 
the King profile and the assumption of the hydrostatic equilibrium of cluster gas. The galaxies in 
clusters, however, constitute a very small fraction of the gravitational mass of the entire cluster due 
to the presence of dark matter. Recent high-resolution N-body/hydrodynamical simulations have 
strongly suggested that dark halos of cluster scales are described by a family of fairly universal 
density profiles. Navarro, Prenk & White (1996,1997, hereafter NFW) proposed a profile: 

^'dm('^) = 7 -, ; -, To; (2) 

(r/rs)(l rlrsY 

where /Oco is the critical density of the universe at z = 0, and and are the concentration 
parameter and the scaled radius whose explicit fitting formulae as a function of the halo mass M 
and cosmological parameters are found in Navarro, Frenk & White (1997). 

This enables one to predict the profiles of the gas and the X-ray surface brightness. Our 
previous work (Makino, Sasaki & Suto 1998, hereafter Paper I) examined the case of the NFW 
profile proposed, and found that the resulting gas density profile is very close to the empirical 
model (0). 

Higher-resolution N-body simulations by Fukushige & Makino (1997), however, indicate that 
the inner density profile of the halo is much steeper than the NFW profile above. This conclusion 
was confirmed later by a series of systematic N-body simulations by Moore et al. (1998). Evans 
&; Collett (1997) show that under some constraints the density profile oc r~^l'^ becomes the 
stationary solution to the collisionless Boltzmann equation. Furthermore, in Paper I we assumed 
the isothermality and the possible effect of the temperature profile was neglected. The present 
paper explores these effects in more details and provides several useful working formulae for 
the X-ray surface brightness profiles as well. In what follows we describe various theoretical 
formulation and do not attempt to compare with real observational data. The comparison with 
observations in this context has been performed by Makino & Asano (1998), Tamura (1998), Xu 
et al. (1998) and Markevitch et al. (1998). While all of their results indicate that the current 
observational data are consistent with the predictions in the universal density profile, the current 
data quality in the spatial resolution of X-ray surface brightness and temperature profiles is 
not sufficient in distinguishing from the empirical /3-model in an unambiguous manner. In this 
respect, upcoming X-ray missions including AXAF and XMM should definitely provide us the 
data suitable for that purpose. This is mainly why we do not attempt any tentative comparison 
with the currently available data, but rather present various theoretical predictions which will be 
useful for future data analysis. 
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The plan of the paper is as follows; in §2 we describe a new series of analytical solutions for 
a family of density profiles of dark matter halos generalizing Paper I. Then we examine effects of 
self-gravity of the gas in §3, and non-isothermality by adopting polytropic equation of state in §4. 
Finally §5 is devoted to discussion and conclusions. 



2. Isothermal gas and X-ray surface brightness profiles from a family of dark 

matter halo potentials 



We generalize the NFW profile (|2|), and consider a family of density profiles describing the 
dark matter halo: 

where x = r/r<j is the dimensionless radius in units of the characteristic scale r^. Then the total 
mass of dark matter halo within the radius r is given by 



M(r) = ATT5cpcorlm{r/rs), (4) 

with 

mix) = I -^du. (5) 

^ ^ io (1 + u^Y 

If one neglects the gas and galaxy contributions to the gravitational mass, the gas density 
profile Pg{r) in hydrostatic equilibrium with the above dark matter potential satisfies 

fcTg dlnpg _ GM{r) 
PgiTLp dr 

where pg and nip denote the mean molecular weight of the gas (we adopt 0.59 below) and 
the proton mass, and we assume that the gas temperature Tg is constant over the cluster 
(non-isothermal cases are considered in §4). Equation (|6|) can be formally integrated to yield 

ln^ = -i?r^"!^dx, (7) 

PgO Jo X'^ 

where 

^ ^ A-KGpgmp5cPcorl 

Note that B is 276/2 in terms of the parameter b defined in Paper I. Also B is rewritten as 

7m(l) Tg ^ ^ 

in terms of the virial temperature defined as 

T^Ur) - ,^^^£|^ (10) 
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with 7(ss 1 — 2) being a fudge factor (see Paper I). 

With the density profile of the form (^), equation converges for /x < 2, and is rewritten as 

Pgi^) = Pgoew[-Bf{r/rs)], (11) 

where 

' Jo Jo {l+u'')^ xJo (l + u^)^ ^ ^ 



In some specific cases, equations (|5|) and (|12D are analytically integrated: 

(i) ^ = 1, = 1, A = 2 (NFW) 

m(x) = ln(l + x) —, (13) 

1 + X 

f{x) = l--ln(l+2;). (14) 

X 

(ii) n = 3/2, = 3/2, A = 1 

m{x) = ^ln(l+x2/^), (15) 
fix) = ^ln(l + ^/^)+ln(l + ^) + -|tan-i^^^ + ^vr. (16) 

(iii) /i = 3/2, 1/ = 1, A = 3/2 



m(x} = 2\n{^Jx + VTTx)-2J- , (17) 

" 1 + a; 



fix) = 2ji^--ln(V^+^/rT^). (18) 
V X X 

In what follows, we focus on the case with /x = a, = 1, A = 3 — a (1 < a < 2). Then 



equation (12) reduces to 



fix) = / 7 — du / -^5 — du. (19) 

■'^ ^ Jo (l + n)3-" xJo (l + n)3-° ^ ^ 

If we set a = 1, equation ( p!9[ ) reduces to the case (i) which corresponds to the original proposal by 
Navarro, Frenk & White (1996,1997) and is worked out in Paper I. Since the power-law slope of 
the inner region is very sensitive to the mass resolution limit of numerical simulations (Fukushige 
& Makino 1997; Moore et al. 1998; Melott et al. 1997; Splinter et al. 1998), we explore the similar 
profiles by changing a. Since it is unlikely that the mass resolution limit affects the asymptotic 
outer halo profile oc r~^, we choose A = 3 — a so as to reproduce the asymptotic behavior. 
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2.1. gas density profile 



For 1 < a < 2, we numerically integrate equation ( ]19| ) to compute the gas density profile 
which is proportional to [-F(x)]^ where we define 



F{x) = exp[-/(x)]. 



(20) 



Figure |i| plots the F{x) for a = 1.0, 1.4 and 1.6 together with the empirical fit to the following 
function: 

1+ 



Xc{a) 



For X <C 1, equations ( p^ to ( [2l| ) are consistent if 

q{a) = 2 — a 

and 



(21) 



(22) 



X- " = {3 - a){2 - a)p. (23) 

We adopt the relation (^2|), but still keep p and Xc as two independent parameters to fit the F{x) 
in the range of 0.05 < x < 5. The results are plotted in solid lines in Figure |l|. For 1 < a < 1.8, 
we find the following empirical fitting formulae: 



Xc{a) = 0.015(2 - a)-2-5 + 0.47(2 - a)°•^ 
as plotted in Figure |2[ 



p{a) = 0.33(2 -a)-l•^^ 



(24) 



2.2. X-ray surface brightness profile 



From the observational point of view, it is more useful to compute the X-ray surface brightness 
profile on the sky: 

1 r°° I 

^ (TT^ (^^iT^Hi^OH\{z) + /2) dl, (25) 

where z is the redshift of the cluster considered, ax is the X-ray (either bolometric or band-limited) 
emissivity, n^{r) is the electron number density corresponding to the profile (pd|), dA{z) is the 
angular diameter distance, and the integration is performed over the line-of-sight at an angular 
separation 6 from the center of the cluster. 

Substituting the isothermal density profiles ([T9|) and (|20|) , equation (p5[) reduces to 

sx(^) = ^"""ffft sie/e,), (26) 



where 6s = rs/dA{z), and 



SicP) . / ^-^^^^E^dx. (27) 
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Throughout the present paper we consider the X-ray bolometric emissivity, and truncate the 
gas density at r = 20rs, i.e., ng{r) = at r > 20rs in evaluating the surface brightness to avoid 
unphysical divergence. This is because 20rs is typically larger than the virial radius and the gas 
profile should not be extended at the larger scale. We integrate equation ( pTj ) numerically for 
{a,B) = (1.0,5.0), (1.0, 10.0), (1.6, 10.0), and (1.6, 20.0), and the results are plotted in Figure |. 
They cannot be reasonably fitted to the conventional /3-model: 

where f3 is the slope parameter, (j)c,i3 is the core radius, and S{0) is the central surface brightness. 
Due to the strong concentration of the dark matter halo in the present models, the resulting gas 
density profiles are also more concentrated than the conventional /3-model. In fact, the deviation 
from the /3-model becomes more appreciable for a larger a. Rather we find that the following 
generalized shape: 

with 

cl)c = 0.3 - l) > (30) 
e = 0.41- 5.4(2 -a)6 + (0.585 + 6.47a-^-^^)B-"'/30^ (3^) 
rj = -0.68 - 5.09(q-1)^ + (0.202 + 0.0206a'^)S^-\ (32) 
provides an excellent fit for 5 < i? < 20 and 1.0 < a < 1.6 in the range of 10~^ < </> < '/'max where 

5(0max) = 10-^5(0). 

Provided that the assumption of the isothermal gas is valid for actual clusters, our fitting 
formulae above can be used directly to probe the shape of the dark matter halo; fitting the 
observed X-ray surface brightness to equations (29) to (32) would result in the values of q, 



and B from the three fitted parameters ^s^o C find r/. In addition if the gas temperature Tg is 
measured from X-ray spectroscopic observations, the concentration parameter 6c is determined via 
equation (^). 



Incidentally let us note that our fitting formulae ( 2£ ) are not unique in the sense that the three 



parameters (pc, C and i] are correlated and written in terms of the two independent parameters 
a and B. In particular, (j)^ and ^ are strongly correlated, and we would obtain an equally good 
fitting formula by changing these two appropriately. On the other hand, r/ is relatively insensitive 
to the choice of 0c or ^. 

Applying the procedure described above to three lensing clusters (A2163, A2218 and RX 
J1347. 5-1145), Makino & Asano (1998) showed that the dark matter profile with a < 1.4 is 
consistent with the ROSAT HUl X-ray surface brightness profiles, and that a ~ 1.4 is preferred 
in order for the X-ray mass estimate to be consistent with their giant arcs. Tamura (1998) also 
reached the similar conclusion on the basis of ASCA/ROSAT observation of A1060. 
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3. Effect of self-gravity of the gas density distribution 

When one properly includes the self-gravity of the gas distribution, equation reads 
kTg dlnpg _ GMtot(r) _ AttG 



fignip dr 



u [pdm{u) + Pgiu)] du. (33) 



With the profile (^) for dark matter halo, the above equation is rewritten in a non-dimensional 
form as 

^ = -4 [—-1—^ + Re^^^A du, (34) 

where we introduce 

g{x) = ln[pg{x)/pgo] and R = Pgo/{ScPco)- (35) 
For p = a, 1^ = 1 and A = 3 — a, the lowest-order perturbation solution for equation (^) at x <C 1 
is o 

In Figure 0, we consider the effect for the NFW profile (a = 1) only, and numerically integrate 
the second-order differential equation: 



d^gjx) ^ 2 dgjx) ^ ^ 
dx"^ X dx 



x{l + xy 



+ 



0, (37) 



where we rewrite equation ( p4D with the boundary condition specified from equation ([Sq). The 
range of R is chosen so that the resulting gas to dark halo ratio: 

^ Mom(x) ln(l + x) - x/(1 + x) io 

approximately ranges 10~^ 10~^ at x = 10. 

We plot the gas density profiles including the effect of self-gravity in Figure ^a. With 
increasing the gas fraction or R, the gas density profile becomes steeper at large radii where the 
self-gravity of gas becomes significant compared with that of dark matter and confines the gas 
itself more strongly. Also as B increases, the gas distribution becomes more centrally concentrated 
since the gas temperature Tg and therefore the pressure gradient against the gravity becomes 
smaller for larger B. 

In Figure ^ we present the gas mass fraction (= Mg/[Mg + Mdm]) for various values of 
B and R. At x ^ 1, puM is proportional to while pgas is approximately constant. So the 
gas fractions of clusters in the present models increase roughly in proportional to x in the inner 
region. For larger x, on the other hand, the behavior is sensitive to the values of B and R. The 
observed baryon fraction of 0.1 ~ 0.2 of clusters indicate that {B,R) = (10,5) and (5, 1) fall in 
an observationally relevant range. In this case, the effect of gas self-gravity is not significant, but 



- 8 - 



cannot be fully neglected either. Also note that the gas mass fraction does not level off anywhere. 
Since the gas mass fraction provides an important constraint on the density parameter (e-g-, 
White et al. 1993), this indicates the need for the quantitative comparison on the basis of the 
numerical integration of the above equation. As discussed in §6, this methodology is feasible under 
the generalized halo potential model. 

4. Example of non-isothermal distribution: polytropic equation of state 

In the discussion above, we have assumed the isothermal gas distribution. Although this is 
regarded as a reasonable approximation to the actual clusters, it is true that some clusters do show 
non-isothermal features. In particular, Markevitch et al. (1997) reported that the temperature 
profiles of clusters appear remarkably similar, and are approximately described by a polytrope 
with the polytropic index T = 1.2 — 1.3 assuming that the gas density profile is given by the 
/3-model. In this section, we consider the effect of non-isothermal gas distribution on the basis of 
the polytropic equation of state for definiteness. More specifically, we adopt the following form for 
the gas pressure P: 



where the subscript denotes the value at the center (x = 0), and the polytropic index T is equal 
to 1 + 1/n. 

As before we assume that the gas is in hydrostatic equilibrium: 




(39) 



or equivalently for the gas temperature: 




(40) 



1 dP _ GM{r) 
Pg dr ' 

and neglect the self-gravity of the gas density. Define the function e(x): 



(41) 




(42) 



then equation (p]) is written in a dimensionless form as 




(43) 



where the constant is defined as 

^ n + \ kTgQ 
Note that as in equation (^) Bp is related to Tvir(r's) as 



n + 1 



B 



(44) 




(45) 
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Equation p^ ) is integrated using the function f{x) (eq.[^) to yield 

e{x) = l-Bpf{x). (46) 

Throughout the rest of this section we consider the NFW profile only just for illustration, but the 
discussion below can be easily generalized to other profiles. 

The temperature and gas profiles in the present polytropic model are determined by specifying 
the two parameters, -Bp and n. It should be noted that the upper limit on Bp is set by the 
maximal extension of the cluster gas; if the gas extends up to a; = Xmaxj then 

Bp < -Bp,i„ax = l//(a;max). (47) 

In the case of the NFW profile, i?p,max are 1.32, 1.18 and 1 for Xmax = 10, 20 and oo, respectively. 
With equation (|45|), the above condition ( p7| ) can be translated to the lower limit on the central 
gas temperature as 

(n + l)7i?p,^axm(l) ~ (n + l)7 ■ ^ ' 

The resulting temperature and gas profiles are plotted in Figure |5|. Note that Tg{x)/TgQ is 
determined by Bp and the halo shape parameters only and is independent of the value of n, while 
the resulting gas density profile is sensitive to n. The gas temperature in the present models starts 
to decrease appreciably around at x ^ 1 (solid lines in Fig.||a), where the density profile departs 
significantly from that in the corresponding isothermal model (dashed lines in Fig.|5|b). The virial 
temperature of the halo T^rix) is significantly different from the gas temperature Tg(x) in the 
present models (also in the isothermal model) indicating the importance of the detailed numerical 
simulations (e.g., Eke, Cole & Frenk 1998; Bryan & Norman 1998; Yoshikawa, Itoh, & Suto 1998) 
in predicting the temperature profile of clusters even in the present context. 

Figure ^ shows the mass profile and gas mass fraction for polytropic models in comparison 
with those for isothermal cases. The gravity of gas mass is neglected, and for definiteness we 
choose the gas mass to be 10 percent of that of the dark halo at a; = 10 which roughly corresponds 
to the virial radius in the NFW model. 

The X-ray surface brightness and the emission-weighted temperature profiles on the sky, 
Sx(^) and Tx{9), are defined by equation (p5|) and 



Tx{9) ^ ^ , , (49) 

/ ax{Tg)nl{^e^dl{z) + P) dl 

J —oo 

respectively. If we consider the bolometric thermal bremsstrahlung emissivity only, the polytropic 
models which we described above yield 

X [e(x)]2"+i/2 

Sx(0) _ JcP ^/X^ - 02 ^^^^ 



^^(0) r [6(X)]2"+V2 ax 

Jo 
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Tx{<P) _ J4> V^^^ Jo ^^^^ 



rx(o) .00 . ieix)r^y^ r ^^^^^^2n^s/2 ^''^ 

where = 9/6^ (§2.2). These profiles are plotted in Figure]^. The dotted curves in the upper panel 
indicate the best-fit /3-model to the surface brightness. It is somewhat surprising that the X-ray 
surface brightness profiles for some models, e.g., with {Bp,n) = (1.0,7.0), are well approximated 
by the /3-model despite the fact that the cluster is far from isothermal. It should be also noted 
here that the projected emission- weighted temperature Tx{9) systematically differs from the gas 
temperature Tg(r) evaluated at 6 = r /dx (e.g., Figs.^ and 0b). 



5. Discussion and conclusions 

In this paper, we have presented a physical methodology to confront the dark matter halo 
mass distribution with the observed X-ray surface brightness profiles of the galaxy cluster. Unlike 
the previous phenomenological prescriptions including the isothermal /3-model, this approach 
enables one to determine the dark matter halo profile directly from the observational data. 
This works in an especially straightforward manner when the cluster is well approximated as 



isothermal; then our fitting formulae (eqs.|24|, to |3^) explicitly links the gas density and 
X-ray surface brightness profiles to the underlying dark matter halo potential as long as the halo 
is described by a family of profiles (eq.^). In §4, we have described a prescription of computing 
the gas, temperature, and X-ray surface brightness profiles for clusters with polytropic equation 
of state. We confirm that the self-gravity of the gas does not affect the density profile significantly 
(§3) in the halo profiles considered throughout the paper. Incidentally Bardelli et al. (1996) 
claimed that the gas mass fraction in A3558 (Shapley 8) is ~ 0.7 at the Abell radius (assuming 
Hq = 50kms~^Mpc~^). If confirmed, the X-ray surface brightness profile for such clusters should 
be computed with properly taking account of the self-gravity of gas, which is in fact feasible as 
shown below. 

When both the self-gravity of gas density and the polytropic equation of state are taken into 
account, our procedure described in the present paper is generalized as follows; 

(i) select a dark matter halo density profile parameterized by (/i, i/. A): 

(ii) fix Bp (eq. ^^), R (eq.|^) and the polytropic exponent n. 

(iii) solve the following equation for e{x) 

0, (53) 



dx"^ X dx ^ 



x^(l + x 
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with the boundary condition at x <C 1: 

(iv) then one obtains Tg{x) = Tgo e(x) and Pg{x) = pgQ [e(x)]", and finally can compute the 
corresponding X-ray surface brightness T,x{0) which should be compared with an appropriate 
sample of X-ray clusters. 

(v) repeat the above procedure to determine the set of parameters /j,, A, n, Bp, and R. The 
latter two quantities are combined to yield the amplitude of the halo density 5cPcO- 

With the current and future spatially resolved X-ray surface brightness and temperature 
profiles of several clusters of galaxies with Einstein, ROSAT, ASCA, AXAF, and XMM, the 
present methodology will be a useful tool in revealing the shape of the underlying dark matter 
halo potential which one has not been able to derive, and thus provide important information on 
theories of cosmic structure formation. 
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Fig. 1. — Gas density profiles predicted from a family of dark matter halo profiles. Solid triangles, 
open circles, crosses, and open squares indicate the results of numerical integrations for a = 1.0, 1.4, 
1.6, and 1.8, respectively. Solid lines represent the best-fits to equation (|2l|) with using equation 
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Fig. 2. — The best-fit parameters of Xc and p as a function of a. Solid lines represent the fitting 
formula (eq.[^4|) which is accurate for 1.0 < a < 1.6. 
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Fig. 3. — X-ray surface brightness profiles predicted from a family of dark matter halo profiles. 
Filled triangles, open circles, crosses, and open squares indicate the results of numerical integrations 
for {a,B) = (1.0,10.0), (1.0, 20.0), (1.6, 5.0), and (1.6, 10.0). Solid lines represent the best-fits to 
equation (p9|) with equations ( pO] ) to (p2D, while dotted lines indicate the best-fits to the conventional 
/3-model (eq. [p^] ) . 
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10-2 lQ-1 1 IQl 

x = r/ 



Fig. 4. — Effect of self-gravity of gas density distribution; a) gas density profile normalized in units 
of the central value; b) gas mass fraction. All models assume a = 1 (NFW profile). 
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Fig. 5. — Profiles of temperatures and densities for polytropic models embedded in the NFW dark 
matter halo potential {a = 1). a) gas temperatures Tg are plotted in solid lines for Bp = 0.5, 
1.0 and 1.5. For reference, the virial temperature Tyir is plotted in dashed and dotted lines for 
{Bp,n,j) = (1.0,7,1.2) and (1.0, 3/2, 1.2), respectively. All values are in units of the central gas 
temperature Tgo; b) gas densities pg for polytropic models (solid lines) and for the corresponding 
isothermal models (dashed lines) with B = {n+ l)Bp in units of the central value pgo- 
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Fig. 6. — Mass profiles for polytropic models embedded in the NFW dark matter halo potential 
{a = 1) with {Bp, n) = (1.0, 7.0), (0.5, 7.0), and (1.0,1.5). a) gas mass profiles for polytropic models 
(solid lines) and for the corresponding isothermal models (dotted lines) with B = {n + l)Bp. The 
gas mass is normalized to be 10 percent of the mass of dark matter halo at x = 10. The halo mass 
profile is plotted in dashed line for reference; b) gas to dark halo mass ratio for polytropic models 
(solid lines) and for the corresponding isothermal models (dotted lines) with B = {n + l)Bp. 
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Fig. 7. — X-ray surface brightness and and emission-weighted projected temperature profiles for 
polytropic models embedded in the NFW dark matter halo potential (a = 1) with (Bp,n) = 
(1.0,7.0), (0.5,7.0), and (1.0,1.5). a) bolometric X-ray surface brightness (thermal bremsstrahlung 
only) profiles for polytropic models (filled circles) and for the corresponding isothermal models (open 
triangles) with B = {n + l)Bp. The dotted lines represent the best-fits to the /? model, b) emission- 
weighted projected temperature profiles for polytropic models (solid lines), gas temperature with 
Bp = 1/2 (dotted line) and Bp = 1 (dashed line). 



